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The localization behavior of the Anderson model with anisotropic hopping integral t for weakly 
coupled planes and weakly coupled chains is investigated both numerically with the transfer matrix 
method and analytically within the self-consistent theory of localization. It is found that the mobil- 
ity edge is independent of the propagating direction. However, the correlation (localization) length 
in the extended (localized) side of the transition can be very different for the two directions. The 
critical disorder Wc is found to vary from t* for weakly coupled planes to t2 for weakly coupled 
chains. 
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The problem of Anderson localization in anisotropic 
systems has attracted considerable attentioru^lj in re- 
cent years, largely due to the fact that the high Tc su- 
perconductors are highly anisotropic. Transport in the 
normal state is metallic in the, ab-plane but appears 
semiconductor-like in the c-axig3. The nature of the c- 
axis transport in high Tc materials is still controversial 
and its understanding may have important consequences 
for the theories of the normal and the superconduct- 
ing statfiL. This paradoxical property has prompted the 
proposalEI that a high T^, material in the normal state 
is actually an insulator, appearing metallic only because 
the inelastic length in the plane is lesajihan the localiza- 
tion length. It has also been arguednQ that a negative 
dp/dT in the c direction alone may signify anisotropic 
localization, with a metal-insulator transition depending 
on the propagating direction, in direct contradiction to 
the predictions of the scaling theory of localizationB. A 
recent diagrammatic calculatioru lent support to such a 
claim. Previous diagrammatic analysisO, however, led to 
the conclusion that the scaling property in anisotropic 
systems remains the same as that of the isotropic sys- 
tems with a simple substitution of the conductance by 
its geometric mean. Given the perturbative nature of 
all the previous work, it is important to carefully study 
the localization behavior of disordered anisotropic sys- 
tems with reliable numerical techniques and to determine 
whetheij-.tlie scaling theory is valid for highly anisotropic 
systemsETn. 

In the present work, we have systematically studied the 
localization properties of a three-dimensional disordered 
anisotropic system described by a tight-binding Hamilto- 
nian with random on-site energies and anisotropic hop- 
ping integrals. We calculate numerically both the con- 
ductance and the localization length of such systems with 
the transfer-matrix methodO. By doing finite-size scal- 
ing analysis on sufficiently large systems, we find the 
metal-insulator transition occurs at the same critical dis- 
order Wc in all directions. Remarkably, at criticality. 



the geometric mean of the ratio of the finite-size local- 
ization length to the width of the bar is independent of 
the anisotropy. This underlies the conformal invariance 
property at the Anderson transition. For systems with 
extremely weak plane-coupling, the critical disorder Wc 
is found to vary with anisotropy t as t*, in good agree- 
ment with a self-consistent theory of localization. How- 
ever, the correlation lengths in the extended regime, as 
well as the localization lengths in the localized regime, 
differ tremendously in different directions. As we will 
discuss later, this difference of the correlation lengths in 
the propagating directions may explain the c-axis trans- 
port of the high-Tc materials. The localization behavior 
of weakly coupled disordered chains, which simulate the 
lD-to-3D behavior, is also studied. The critical disorder 
Wc is found to behave as 1 2 . 

We study the following Hamiltonian for an anisotropic 
3D disordered model 



H = ^e«|n >< 7l| -I- ^ ]tnm 



\n >< m| 



(1) 



where the sites {n} form a regular cubic lattice. The hop- 
ping integrals are nonzero between nearest-neighbor sites 
only and depend on directions, in general, ^ ty ^ tz- 
We normalize all energies by the largest hopping integral. 
For systems with weak plane coupling, the hopping in- 
tegrals are given by {1,1, t}. As a convention, we have 
assigned the direction with the large and small hopping 
integral as the parallel (||) and the perpendicular (_L) di- 
rection, respectively. Disorder is introduced by choosing 
randomly site energies e„ within [-W/2, W/2]. 

The metal or insulating nature of the Anderson Hamil- 
tonian, Eq. (1), can be determined by investigating the 
scaling properties of either the conductance or the lo- 
calization length of finitCj-ipstems calculated from the 
transfer-matrix techniquea3'tll. The conductance G of an 
M X M X M cubfiJs calculated from the multi-channel 
Landauer formulate 
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G{M) = ^Tr{T+T) (2) 

where T is the transmission matrix. For the localization 
length calculation, one considers a bar of length N and 
cross-section M x M. One determines the largest local- 
ization length Am as N ^ cxd from the smallest Lyapunov 
coefficient of the product of the random transfer matrices. 
From a plot of Am/M (or G) vs. M, one can determine 
the localization properties of the system. For localized 
states (W > Wc) the ratio Am/M (or G) is expected to 
fall with increasing M, while for extended states (W < 
Wc) Am/M (or G) should rise instead. At the mobihty 
edge trajectory (W=Wc), Am/M (or G) is independent 
of M and this behavior defines the Anderson transition 
point and the critical disorder W = Wc. In our calcu- 
lations, we have used systems with M = 5 - 17 and N 
was at least 5000. For the conductance calculations, M 
was up to 22. Due to the nonself-averaging nature of 
finite-size systems, an average over many random config- 
urations (up to 500 for the M = 20 case) must be taken 
to suppress the large fluctuations. 

In Fig. 1, we present our numerical results for the 
dimensionless conductance g ~ -^jj^ and the scaled lo- 
calization length Am/M versus M at E=0 for the case of 
weakly coupled planes with coupling constant t = 0.1, 
for both the perpendicular (lower panels) and parallel 
(upper panels) directions. These results strongly suggest 
that for sufficiently large M, the critical disorder defined 
by an A^I-independent g and Am /M converges to the same 
value, Wc ~8.4, for both propagating directions. Notice 
that if only sizes M< 11 were used, which are appropri- 
ate in the isotropic case, one would then have erroneously 
concluded that wl > W^. Another important point is 
that the value of Aj\//M at W = Wc, called Ac, is differ- 
ent in the two propagating directions. We find that A^ 

= tAl for all the t's we have examined. The critical con- 
ductance in the two directions at t=0.1 is approximately 
g^ 2± lO^^gc. This is a much stronger variation thanJthe 
g-^ ~ t^ gll predicted by the diagrammatic analysisB of 
the anisotropic model. This may be an indication of the 
existence of high order corrections to the conductivity 
that is not simply proportional to the bare conductivity. 
This point needs further study. 

In order to extrapolate to infinite system size (M— > oo), 
it is necessary to investigate the scaling behavior of 
Am/M. In the center of the band (E = 0) and close 
enough to the transition point, it turned out to be pos- 
sible to establish a scaling function for both propagation 
directions within the accuracy of our numerical results. 
The scaling function f(x) behaves as 1/x in the x ^ 
limit for extended states (while for localized states f(x) 
~ x). For X ^ oo, f(x) approaches a constant value that 
depends on the propagating direction. For t ~ 0.1, we 

obtained Ac = 1.2 and A^ = 0.12. This suggests that 

A^ = tAl for general t. We have, indeed, checked that 
this formula is correct for any t. Another important rela- 



tion we were able to obtain is the geometric mean of the 
different critical values Am /M in the different directions 
is independent of the anisotropy t. We derived 

(A^A^A^)i/3 ^ 0.6 (3) 

In the case of weakly coupled planes, Eq. (3) becomes 
((Ac)^A^)^/'^ ~ 0.6. This relation may have impor- 
tant consequences for the existence of conformal invari- 
ance at pthe. critical point of the Anderson localization 
problemtjo. 

In Fig. 2, we plot the localization length Lc and the 
correlation length ^ as a function of W for E = and t 
= 0.1 for the parallel and the perpendicular directions. 
Both Lc and f diverge at the critical disorder Wc. How- 
ever, both ^ and Lc differ substantially for the different 
propagating directions. In particular, in the extended 
regime where W<Wc, > For example, at W=5, 
^^=100, and ^" ~1.5, in units of the lattice constant, 
which has been taken to be one. In the localized regime 
where W>Wc, lI > L^. For example at W=10, l!!=55, 
and L|f ~5. Our numerical results approximately fol- 
low the theoretical predictioni = t^^^ and = thi. 
However, the exponents are expected to be the same and 
this has been confirmed from our calculations. We find 
that = 1.3 ± 0.1 and z/" = 1.3 ± 0.3, in agreement 
with each other and with v = 1.3 ± 0.1 for the isotropic 
system within the numerical accuracy. 

The difference between and ^""^ is very important 
and can nossibly explain the normal state transport 
properties^ of the high-Tc materials. The correlation 
length ^ measures the strength of the fluctuations of the 
wave functions in the extended regime. For length scales 
larger than ^, the wave function looks uniform, while for 
length scales smaller than the wave function has strong 
fluctuations. The relevant length scale is the inelastic 
mean free path which behaves as T^^*, with probably 
p=l/2. When tin < ^, a phenomenon called incipient 
localization takes place and conductivity |iS| controlled by 
lin- A convenient interpolation formulae^, valid in the 
conducting regime is given by 



where a and b are constants of order unityE£l. For high 
enough temperatures the conductivity is given by the 
regular metallic behavior, where Gph = ajpT/Air and the 
mean free time t T~^. Because Gph and tx/oc describe 
independent physical processes, we can add the corre- 
sponding resistivities. Therefore p — pioc + Pph, where 
p — 1/a. The experimental behavior of the high-Tc ma- 
terials can be understood if we assunae that the high- 
Tc oxides, instead of being insulatorslJ, are disordered 
anisotropic metals with a large anisotropic correlation 
length. At low T, the resistivity is dominated by pioc- 
Once £in becomes shorter than ^ in the perpendicular di- 
rection (c-direction) , Eq. (4) suggests that there is 
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an sharp downturn in the perpendicular resistivity as the 
temperature increases. This trend will eventually stop 
at some T when the regular resistivity begin to domi- 
nate. Then will start increasing linearly with T, as 
in the regular metallic behavior. In the parallel direc- 
tion, ^11 is always smaller than the inelastic length and 
the transport in the plane remains metallic. The above 
analysis neglects correlations and dynamic disorder that 
could also affect substantially the transport propertied. 

The dependence of Wc on the anisotropy t for the weak 
plane coupling has been also numerically calculated by 
the transfer matrix technique. The numerical values of 
Wc can be fitted well with a single power law dependence 
Wc=15.4 t3, for t < 0.8. This is clearly seen in Fig. 3. 
This dependence is in marked contrast to whatj-one would 
expect based on reasonable heuristic argumentlllajid also 
predicted previously by more elaborate theoriesLlij which 
give a much weaker t dependence, Wc ~ 1/ y/\£nt\. The 
t dependence is also different-irom the results obtained 
for the weakly coupled chainsc3, which give Wc '^t2. To 
understand the t* dependence of Wc, we examine the 
localization transiticm starting from the results of the di- 
agrammatic analysia3„ The maximally crossed diagrams 
produced a correctionQ to the zero-temperature configu- 
rationally average conductivity (i=x,y,z) of the form 
given in Eq. (2.3a) of Ref. 4. An equivalent criterion 
in the tight-binding representation for the mobility edge 
is given by Eq. (2.5) of Ref. 4. For the anisotropic 
case we made the choice that the effective lattice con- 
stant tti is proportional to anisotropic mean free path ii. 

£i ~ (< vf >)i/2r = [ < t; >< ^ > Y^^T, i.e., we as- 
sume that T is isotropic. Since aio/af is independent of 
i, the localization criterion is satisfied simultaneously in 
all the directions and can be written as 
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(27r)3 



dq 



1 



Ej=i(2 - 2cosgj) 



(5) 



where n = OLi and 2G|"£,(E=6)=0.505462 is the 
3D Green's function of the isotropic system at the band 
edgalj. Notice the same equation applies to the isotropic 
case except now both the conductivity and the mean free 
path enter as geometric means. 

In the limit of weak coupling t^O, Eq. (5) yields 
T t^i for weakly |eoupled planes and r ~ for 
weakly coupled chainsE2l, and therefore Wc ( taking into 
account that t ~ W^^ ) is proportional to t* and ta , re- 
spectively, in agreement with the transfer-matrix method 
results. We have systematically calculatedc2l Wc vs t, for 
general t, within the CPA, by using Eq. (5). There is 
quantitative agreement between the CPA and the trans- 
fer matrix results. We feel the choice of the effective 
lattice constant a^ being proportional to an anisotropic 
mean free path ii is the proper one and the good agree- 
ment with our numerical results further supports this. 

In summary, we have numerically studied, by the 



transfer-matrix method technique, the localization prop- 
erties of weakly coupled planes. We found that only one 
mobility edge exists for both propagation directions, i.e., 
the states in the planar direction become localized with 
exactly the same amount of disorder as the states in the 
perpendicular direction. However, the correlation length 
^ in the extended side of the transition, as well as the 
localization length in the localized side, can be very dif- 
ferent for the two propagating directions. This behavior 
of ^ can possibly explain the transport properties of high- 
Tc materials. The critical value of disorder Wc is propor- 
tional to 1 4 for weakly coupled planes and is proportional 
to 1 2 for weakly coupled chains. These results are found 
to be in satisfactory agreement with the predictions of 
the self-consistent theory of localization which incorpo- 
rates the idea of length scale rescaling. However, the 
conductance in different directions does not satisfy the 
relation predicted by the diagrammatic analysisB. Work 
to further understand the scaling properties of the con- 
ductance are in progress. 
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FIG. 1. The conductance g (a and b) and Am/M (c and d) 
plotted as a function of M for E=0, t=0.1, and various values 

of disorder W, for both propagating directions. The mobility 
edge is at Wc ~8.4 for both propagating directions. 



FIG. 2. Localization length Lc or correlation length ^, ob- 
tained by finite size scaling, as a function of disorder W for 
E=0 and t=0.1, for the parallel and perpendicular propagat- 
ing directions. Notice that for W<Wc ~8.4, > ^H, while 
for W>Wc, lI! > L^. 



FIG. 3. The critical disorder Wc for obtaining localized 
states at E=0 vs the anisotropy t. 
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